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ABSTRACT 

We explore the realm of post-Asymptotic Giant Branch (post-AGB) stars from a 
theoretical viewpoint, by constructing synthetic population of transition objects, proto- 
Planetary Nebulae, Planetary Nebulae Nuclei, and post-Planetary Nebulae objects. We 
use the Montecarlo procedure to filter out the populations accordingly to a given set of 
assumptions. We explore the parameter space by studying the effects of the Initial Mass 
Function (IMF), the Initial Mass-Final Mass Relation (IMFMR), the transition time 
(ttr ), the envelope mass at the end of the envelope ejection (M^), the planetary nebula 
lifetime tpN , the hydrogen- and helium-burning phases of the central stars. The results 
are discussed on the basis of the HR diagram distributions, on the My — t plane, and 
with mass histograms. We found that: (1) the dependence of the synthetic populations 
on the assumed IMF and IMFMR is generally mild; (2) the M^ indetermination 
produces very high indeterminations on the ttr and thus on the resulting post-AGB 
populations; (3) the synthetic models give a test check for the ratio of He- to H-burning 
PNNi. In this paper, disk post-AGB populations are considered. Future applications 
will include Magellanic Clouds PNe, and populations of bulges and elliptical galaxies. 

Subject headings: Stellar Evolution; Stellar Populations; post-AGB stars; Planetary 
Nebulae Nuclei; Initial Mass Function; Galaxy 

1. Introduction 

The post-asymptotic giant branch (post-AGB) phase of evolution, when stars leave the red 
giant region to eventually become white dwarf (WD) remnants, is characterized by a series of 
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intriguing events and diversities. The beginning of the post-AGB phase is conventionally set at the 
time of cessation of the high mass loss rate episode M > 10~^ M0yr~^ which has almost completely 
stripped the star of its hydrogen-rich envelope. After this superwind (SW) quenching, the central 
star shrinks in radius at almost constant luminosity, thus getting hotter and hotter. Initially it is 
highly obscured by the fossil superwind around it, but in ~ 100 — 1000 yr this circumstellar material 
becomes optically thin and the star returns to be observable also at optical wavelengths. The so- 
called proto-planetary nebula phase (proto-PN) begins at SW quenching and ceases when the star 
becomes hot enough (Tefj > 10, 000 K) to start photo-ionizing the materials that were ejected at 
relatively low velocity (~ 10 — 20 km s^^) during the superwind phase. As cascade recombination 
and free- free emission make the nebular material shining in the optical, the object has turned 
into a planetary nebula (PN). The hot star is now emanating a very fast (several 1000 km s~^), 
radiatively-driven wind which dynamically interacts with the fossil superwind, occasionally shaping 
it in extravagant forms. As nebular expansion continues, the luminosity and surface brightness of 
the nebula drop as a result of its decreasing emission measure and increasing transparency to the 
ionizing photons. A time inevitably comes that its surface brightness becomes too dim for the 
nebulosity to be noticed by terrestrial astronomers, and the still hot post-AGB star is now slowly 
evolving towards its final WD configuration. 

Occasionally, however, before becoming a WD a final thermal pulse of the still active helium 
burning shell injects enough energy into the outer layers to cause a dramatic expansion of the post- 
AGB star, and sending it back towards the red giant region for a while. Depending on the details 
of this metamorphosis, this residual hydrogen-rich envelope can be engulfed by the underlying 
convectivc helium shell, with its hydrogen being then diluted and burned, and the energy thus 
released causing further expansion. A hydrogen-deficient carbon star - a class of objects including 
RCrB stars - is thus formed, but even when the ingestion of the hydrogen envelope does not take 
place, modest wind mass loss suffice to remove virtually all this residual envelope, thus exposing 
a bare core with a helium, carbon and oxygen atmosphere emanating a Wolf-Rayet like spectrum. 
Therefore, depending on whether or not post-AGB stars miss to experience such a final thermal 
pulse, stars enter the WD stage with or without a hydrogen dominated atmosphere, a condition 
that certainly has to do with the observed dichotomy of WDs into the DAs and the non-DAs groups 
(respectively with and without hydrogen). Life and times of AGB and post-AGB stars have been 
extensively reviewed by Iben &; Renzini (1983, IR) and Iben (1993, 1995). 

Besides its interest per se, the post-AGB phase in stellar evolution has also a number of in- 
teresting connections to other important astrophysical issues. For example the PN composition 
provides an essential tool to study the previous nucleosynthesis which has taken place during the 
AGB phase and before, thus offering important input for modeling the chemical evolution of galax- 
ies. Moreover, PNs are used as tracers of stellar populations, e.g., for determining the stellar death 
rate of galaxies, and arc being used as standard candles in cxtragalactic distance determinations. 
Finally, after neutron stars, post-AGB stars are the hottest objects a stellar population can gener- 
ate, with the most massive of them coming close to one million K, though for a short while. Their 
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contribution to the UV and soft X-ray part of the spectrum is therefore of potential importance, 
especially in old populations. 

The broad scenario sketched above for the post-AGB evolution has now gathered fairly general 
consensus - at least for most of its aspects. This picture has gradually emerged thanks to an 
extended set of stellar evolutionary calculations, starting with the pioneering work by Paczyiiski 
(1971), and continuing with a series of other relevant studies (which include: Harm &; Schwarzschild 
1975; Schonberner 1979, 1983; Iben et al. 1983; Iben 1984, 1985, 1987; Iben k Tutukov 1984; Iben 
& MacDonald 1985, 1986; Wood & Faulkner 1986; VassiUadis & Wood 1994; and Blocker 1995). 
Still on the theoretical side, several of the evolutionary transformations mentioned above have been 
first sketched in a scattered series of papers (Renzini 1979, 1981a, 1981b, 1982, 1983, 1989, 1990, 
hereafter PI through P7), along with a number of conceptual tools which may be useful for their 
understanding. This paper is now an attempt at presenting these ideas in an orderly and systematic 
form. To do so we have proceeded with the construction of simulated post-AGB populations, as 
this technique offers an opportunity to see in a rather straightforward way what are the observable 
effects of the various theoretical ingredients, and therefore to estimate the relative uncertainties 
when - proceeding in the opposite direction - one attempts to use evolutionary theory to infer 
astrophysical quantities from the observational data. 

The paper is organized as follows: In §2 we describe our procedure for the construction of 
post-AGB synthetic populations, in a way that makes them easily reproducible. We therein also 
discuss the role of the parameters that are at variation in our population synthesis. In §3 we present 
our results: synthetic post-AGB populations are plotted within the HR diagrams, and My-age 
( My-PN radius) distributions; mass distributions are also shown. Our conclusions are summarized 
in §4, where we also draft future applications. In this paper we deal with disk populations with 
solar composition. Future work will extend the parameter space to other populations as well. 

2. Synthetic post-AGB populations: ingredients 

2.1. Interpolation of the evolutionary tracks 

The stellar models used in this paper are interpolation of Vassiliadis and Wood (1994, hereafter 
VW)'s hydrogen burning evolutionary tracks, calculated for post-AGB masses between 0.569 and 
0.9 Mq, and solar composition (from here onward, M indicates the post-AGB stellar mass). To 
extend the baseline of our population synthesis, we extrapolate to lower and higher masses by using, 
as templates, the log T^s -log L tracks for M = 0.546 and 1.2 Mq respectively from Schonberner 
(1983) and Pacziynski (1971). We choose the VW database for its homogeneity and wide mass 
range. The VW data set includes also some lower metallicity tracks, which will be used in a future 
paper for simulating the Magellanic Cloud post-AGB populations. 

In Figure 1 and 2 we plot effective temperature and luminosity versus time, as read directly 



-4- 



from the VW tracks^. These tracks do not support straightforward interpolation, since they cross 
one another in several points. We divide each track into parts roughly representing physical phases, 
and then we interpolate the appropriate normalized functions for each phase. 

The effective temperature-age curves have been interpolated as follows. The six log Tcfj - 
log t tracks have been divided into four parts, which hereafter arc called (a) H-burning phase, (b) 
quenching phase, (c) cooling phase, and (d) white dwarf phase. In Table 1 we list the characteristics 
of the VW models. Column (1) give the phase, column (2) the mass of the evolutionary models, 
column (3) the characteristic time of each phase and mass. Note that each phase starts at the end 
of the previous one, and that the H-buming phase starts at t=0. Table 1 also indicates the values 
of log Tefj and log L at the end of each phase, in columns (4) and (5). Prom here after, the capital 
letters H, Q, C, and W will flag the physical variables at the end respectively of the hydrogen 
burning, the quenching, the cooling, and the white dwarf phases (e.g., tii{0.5Q9 M0)=3.215 x 10^ 
yr). For each phase, which correspond to a well defined age interval, we have studied the best way 
to normalize the log Tgs -log t function to eliminate the crossing of the tracks of contiguous masses. 

In Figure 3 we show the temperature evolution on the hydrogen burning phase.'^ Here 
log Teff has been plotted against the normalized age 

tnn = t/tji; (1) 

the dots represent the normalized evolutionary data points. The interpolation of the normalized 
function is straightforward in this phase. Similarly, Figures 4, 5, and 6 show respectively the 
log Teff curves for the quenching, the cooling, and the white dwarf phases. In the quenching phase 
(Fig. 4), it is possible to interpolate directly log T^g versus tnn- The normalized age in the quenching 
phase is: 

tnn = {t-tn)/{tQ-tH). (2) 

In the cooling phase (Fig. 5), we use the interpolating function: 

Tnn = log Teff /log T,s,Q (3) 

with normalized age: 



^In the VW database, all log T^s -log t evolutionary tracks start at t=0, corresponding to log Teff =10,000 K. The 
curves in Fig. 1 start from the second evolutionary data point for clarity. 

''it is understood that all temperatures are in Kelvins, all masses in Mq, the stellar luminosities in Lq, and the 
ages in years, unless otherwise noted. 
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tnn = {t- tQ)/{tc - tq). (4) 

In the last branch of the evolutionary track, the white dwarf phase (Fig. 6), we normalize the 
effective temperature as: 

rnn = Iogreff - log Teff,c, (5) 

where the normalized age assumes the values of: 

Un = it- tc)/{tw - tc) (6) 

(note that of Eqs. 3 and 5 do not correspond to an actual model temperature). 

Let us explore the luminosity tracks, for which we use a a similar interpolation approach. The 
luminosities against normalized ages are plotted, for phase H, in Figure 7. The luminosity interpo- 
lation of this phase is evidently straightforward. The quenching phase is very well interpolated by 
a straight line, calibrated for each mass as 

Lnn = A + B [{t - iH)/(tQ - tu)], (7) 

with A = 0.041 and B = 1.02. The cooling phase normalized luminosity is plotted in Figure 
8 for tracks corresponding to masses between 0.569 and 0.754 Mq. For larger masses, a better 
interpolation is achieved by using the power function, with exponent 0=0.23 for 0.754 ^ Af < 0.8 
Mq, and a=0.26 for M > 0.8 Mq. Figure 9 shows the white dwarf phase of the luminosity (vs. 
time) curves, where log L has been plotted against the normalized time, 

tnn = it- tc)/{tw - tc). (8) 

Even in this case, the interpolation is straightforward. 

After breaking the log T^s and log L tracks into phases, and constructing the normalized 
functions as described above, our procedure includes the following steps, for each given mass and 
age. (a) We identify the mass-interval to be used. For example, if M = O.8M0, we will use the 
two enclosing tracks for the interpolation, corresponding to 0.754 and 0.9 Mq. (b) We identify the 
phase of the evolution to consider, given the age and mass. To do so, we evaluate in, tq, tc, iw for 
the given mass. Let us suppose, in our example, that t = 400 years. We obtain respectively 775, 
859, and 29,689 yr for the H, Q, and C characteristic times. This means that at t = 400 our 0.8 
Mq star is in the H-burning phase, (c) We calculate the normalized t iniG, ^nn* our exctmple, 
^nn = 0.52. (d) On the appropriate tracks (in our example, Figs. 3 and 7) we read off log T^q and 
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log L in the appropriate range of t and M. In our example, in Fig. 3 we draw the vertical line 
X = 0.52 and read off the log T^s values from the two upper tracks, the four points immediately 
before and after the vertical line, (c) We interpolate, on the higher and lower mass tracks, and find 
the log Teff and log L that intersect the x = tnn vertical line; (4) we interpolate vertically versus 
log M, where M is the given mass. 

When evaluated on the masses of the original VW evolutionary tracks, the synthetic temper- 
ature and luminosity curves are indistinguishable from those plotted in Figures 1 and 2. We have 
studied the interpolation offsets in detail. The relative errors in the linear temperatures, (Tgff gyn- 
Teff,vw)/Tefr,vw ^1 result to be collimated between -0.005 and 0.005 in most evolutionary phases, 
with the exception of the very fast quenching phase, where the linear errors on the temperature 
are rather between -0.02 and 0.02 for the lower 5 mass tracks, and between -0.02 and 0.1 for the 
0.9 track. The relative errors in the luminosities, (Lsyn-Lvw)/LvW) are between -0.01 and 0.01 in 
most evolutionary phases, with the exception of the quenching phase, where the linear errors on 
the luminosity are between -0.15 and 0.15 for the lower 5 mass tracks, and between -0.15 and 0.7 
for the 0.9 track. 

When using the interpolated values au lieu of the evolutionary ones, we obtain a perfect sub- 
stitution. In the quenching phase, the substitution would displace the temperatures by a marginal 
fraction, and the misplacement of the luminosity curves would propagate only in a few percent 
error on log L. In all other cases, the substitution reflects in virtually no errors. We conclude that 
the interpolations are excellent and can be reliably used as evolutionary tracks for all stellar masses 
between 0.569 and 0.9 Mq. 

We obtain extrapolations to higher and lower masses using the log Tgg -log L templates of 
0.546 and 1.2 solar masses from Schonberner (1983) and Paczinski (1971). Only the shape of the 
log Tes -log L relation is used here, not the actual evolution of the physical parameters, which are 
instead extrapolated directly from the 0.596 and 0.9 Mq tracks by VW, for homogeneity. 

2.2. The initial mass function 

The synthetic distribution of post-AGB stars (e.g., on the HR diagram) depends on the adopted 
IMF and history of star formation. Note that a given distribution of post-AGB stars obtained with a 
certain IMF and star formation (SR) history can be reproduced also with a different IMF provided a 
properly tuned SF history is adopted. We have renounced to play with two independent parameters 
- such as the slope of the IMF and the e-folding time of the SF rate - and have exclusively explored 
the effect of changing slope of the actual initial mass distribution 'ip{Mi) oc MC^^'^^\ Therefore 
X is effectively the slope of the IMF only if a constant star formation rate is assumed. This may 
not be a bad approximation for the Galactic Disk (Scalo 1998). In our simulations it is supposed 
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that only stars in the mass range 0.85 < Mi < 9Mq experience the thermal pulses on the AGB 
(TP-AGB) and the post- AGB phases. The upper limit of the mass domain has been chosen in the 
light of the carbon ignition limits described in Iben (1995). We include in our investigation the 
classic Salpeter's (1955) IMF, with constant index for each mass range. Furthermore, we inspect 
the effects of a IMF with variable exponent, such as the one by Miller &; Scalo (1979), the low-mass 
extension by Kroupa et al. (1991), and the IMF recently proposed by Scalo (1998). In all three 
cases, the function index changes for masses smaller or larger than Mi= IM©. 

2.3. The initial mass-final mass relation 

The initial mass-final mass relation (IMFMR) may have a major impact on the resulting syn- 
thetic distributions. To illustrate the case we have constructed simulations adopting four different 
IMFMRs, either from theoretical calculations or from observations: (1) the old theoretical Renzini 
&; Voli (1981) relation as analytically approximated by IR ^; (2) a more recent theoretical IMFMR 
introduced by Ciotti et al. (1991), which is an attempt at incorporating new important results of 
evolutionary calculations; (3) the classic empirical IMFMR proposed by Weidemann (1987) on the 
basis of WD masses observed in open clusters; and (4) a new empirical relation presented by Her- 
wig (1997) that includes new observations of cluster white dwarfs in the Pleiades, the Hyades, and 
NGC 3451. Figure 10 shows the four IMFMRs (note that in this paper we use M=Mf). 

A brief justification for the most recent theoretical relation (Ciotti et al. 1991) is appropriate. 
The old theoretical IMFMR (IR) was obtained interpolating on an insufficient grid of stellar models, 
and assumed a universal core mass- luminosity relation. Lattanzio (1989) pointed out that the 
straight line approximation for the core mass at the first AGB thermal pulse (as a function of 
Mi, and for Mi < SMq) was much steeper than resulting from actual evolutionary calculations. 
Furthermore, Blocker Sz Schonberner (1991) have found a major breakdown of the core mass- 
luminosity relation (a key ingredient in theoretical IMFMRs) for those stellar models in which the 
so-called envelope burning process is activated, i.e., for Mi > 3Mq (cf. Renzini & Voli 1981). In 
practice the former effect results in a flattening the IMFMR below ^ SMq, the second one in a 
flattening of the IMFMR above ~ 3M0. The new theoretical IMFMR incorporates these flndings, 
while keeping the same parameterization of mass loss processes (wind and superwind) as adopted 
by Renzini &; Voli, with = 0.5 and 6=1. 

On the empirical side, the new relation by Herwig (1997) is not very different than the new 
theoretical relation. In addition, we should mention that new data points by Jeffries (1997) of 
white dwarfs in NGC 2516 also lie very close to Hcrwig's relation. Furthermore, it should be noted 
that more recent calculations presented in conferences (e.g., Lattanzio & Forestini 1999; Blocker 
1999) show that the envelope burning occurs for M>4 Mq, bringing the theoretical MIMFR hardly 



In this paper we flag this IMFMR with 'IR'. 
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distinguishable from Herwig's (1997). 



2.4. The transition time and the post-AGB envelope mass 

At the TP-AGB phase, the red giant envelope is ejected by the superwind. The subsequent 

evolution of the star to PNN depends in a substantial way on the amount of envelope mass left on 
the star at the SW quenching, M^. No matter how the transition between the AGB and the PN 
illumination occurs, the remnant envelope mass at the end of the SW plays a central role. This is 
just unfortunate when trying to simulate this evolutionary phase: in fact, Mf- is not defined by 
stellar evolution (given the hydrodynamical nature of the superwind), and only approximations or 
guesses can be made about its entity and its possible relation to the physical parameters. can 
indeed be considered a free parameter, in the sense that there is no theoretical nor observational 
constraint to fit it, nor to give it an exact dependence on any nebular or stellar parameter. 

In principle, the higher the core mass, the higher the stellar luminosity, and the lower the 
envelope mass left on the AGB star. But this also depends on the thermal pulses that occur at 
the TP-AGB, and it is hard to know from first principle how the superwind ejection changes the 
stellar structure. We examine later in this section how we choose Mf- for our simulation. 

To understand the transition between the AGB and the PN phases, we should introduce the 
timescales involved in the different phases. Immediately after the superwind quenching, the mass 
loss due to stellar wind dominates, until the star detaches itself from the AGB. The time scale at 
which this phase occurs is the wind time scale, and can be written as: 

u = (9) 

where AMc is the difference between the residual envelope mass at the superwind quenching, 
M^, and the envelope mass at the detachment from the AGB (i.e., at a later phase), M^. The 
wind mass loss rate (MLR) is taken from Reimers (1975), evaluated at the AGB temperature of 
5000 K. Reimers' approximation seems reasonable in the considered evolutionary phase, that is, 
after the superwind quenching and before the onset of the radiation-driven wind (see Blocker 1995). 
In effect, the MLR right after the detachment from the AGB declines as the stellar temperature 
increases (Blocker 1995); but in the present application, i.e., to evaluate the transition time, this 
choice does not affect the results. 

For larger temperatures, in the blueward runaway, the mass consumption occurs at a nuclear 
timescale, following the Equation: 

«n = ^ ■ (10) 
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AMe is now the difference between the smaller among Mf' and M^, and the amount of 

envelope mass remaining at the illumination of the nebula, M^; the luminosity in the equation is 
the plateau luminosity ^; and the variables and En are the envelope hydrogen abundance and 
the energy released by the nuclear burning of one gram of hydrogen, respectively. 

Finally, the transition time will scale with the shorter among the wind and the nuclear time 
scales. We set the reciprocal of the transition time to be the sum of the reciprocals of the wind 
and the nuclear times. 

The transition occurs on a thermal timescale in the case in which M^< , thus in the case 
in which either the nuclear or the wind time is zero. The formula to be used then is: 



ith = ^^^If^, (11) 
where L and R are the stellar luminosity and radius (see also P6). 

It is worth recalling that the synthetic tracks obtained in §2.1 have their zero age points at 
log Teff =4.0, which are approximately reached when the transition has been completed. At 10,000 
K, the star is able to ionize hydrogen, although the complete nebular transparency at the optical 
wavelengths occurs generally at slightly higher temperatures (Kaufl et al. 1993). In all events, we 
use the synthetic tracks from their natural zero point. 

Evolutionary models for post-AGB stars predict the amount of mass available in the stellar 
envelope as a function of the effective temperature (Schonberner 1983, Paczynski 1971). We use 
Schonberner's models (M< 0.65) to obtain M,? (envelope mass at the detachment) and M|?^(envelope 
mass at nebular illumination, when Teg =10,000 K) as functions of M, by linear interpolation. For 
M>0.65 Mq we obtain and by scaling Paczynski's 1.2 M© model to Schonberner's 

values. 

As discussed at the beginning of this section, Mf- is quite unconstrainted, and its value affords 
unpredictable variations. We set ourselves to explore a wide parameter space with three different 
assumptions for M^. 

First, since the evolutionary models have shown that the envelope mass at several evolutionary 
phases in the post-AGB is inversly proportional to the core mass, we could envision the possibility 
that also at the superwind quenching has an inverse dependence on M. We parametrize this 

case following the usual models as guidelines, by taking as the envelope mass at the onset of 
the post-AGB models. Figure 11 illustrates how the characteristic timescales vary if M^=f(M). 
The transition time peaks for the very small mass models, and then it declines for larger masses. 



'^The plateau luminosity, i.e., the luminosity in the early post-AGB phase, is Lpl = 56694 x (M — 0.5) for VW 
models. 
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Second, we chose a constant Mf-, independent of mass. In Figure 12 and 13 we analyze the 
consequences of assuming two values of M^, M^=3xl0~^ and M|?=5xl0^'^. In the first case 
(Fig. 12), the transition time follows up close the thermal time scale, then the nuclear time scale for 
m>0.85 Mq, to decline for larger masses. In the second case (Fig. 13), the transition time rises to 
almost 6000 yr for M~0.55 Mq. These two cases with constant Mf- are illustrated only for the 
sake of showing extreme cases, but it is very unrealistic that all stars end the unstable SW phase 
with exactly equal Mf-. 

Third, we use random values of for our calculations. In Figure 14 we show the transition, 
wind, and nuclear time scales for random M^, with maximum equal to Mf'=5 x 10"^. The 
transition time peaks at about 0.6 solar masses, then declines. The same simulation with different 
maximum gives similar results, except for the vertical scale. Naturally, giving the randomness 
of the extraction, the population plotted in Figure 14 is one of the infinite possible extractions with 
random lesser than 5 x 10~^ Mq. 



2.5. The final helium-shell flash 



During the post-AGB phase hydrogen is burned in a shell, and therefore the mass of the 
helium buffer zone between the C-0 core and the hydrogen shell keeps increasing. In some cases 
the increase of the buffer mass can lead to a last thermal pulse (also called flash) of the helium 
burning shell (Schonberner 1979). This happens if the star left the AGB with a sufficiently massive 
buffer zone, so that its further increase during the post-AGB allows to reach the critical value for a 
last flash to erupt. This chance is therefore related to the phase in the AGB thermal pulse cycles at 
which the SW envelope ejection takes place. Detailed model calculations show that as a response 
to the final flash stars undergo an extended loop in the HR diagram. The model star would expand 
back to the AGB very rapidly, then it would evolve again as a post-AGB star, powered by helium 
burning; the duration of the He-burning phase is about three times longer than the fading time of 
hydrogen-burning post-AGB stars (Iben 1984). These calculations also show that during the power 
down phase the rate of luminosity decline is nearly constant, as opposed to the case of hydrogen 
burning post-AGB stars in which an abrupt drop follows the plateau phase. 

To incorporate the effect of the final helium-shell fiash (FF) we have proceeded as follows: we 
assume that such stars return to (L, Tgfr) = (-LpL, 10, 000 K) upon a final flash, and their luminosity 
fades then linearly with time to Lp in a time 3 x (where tp is the total fading time for the 
H-burning star, see § 2.7.2), and Lp is the luminosity of the corresponding H-burner, i.e.: 



L = LpL - ~ ■ t'. (12) 

We then locate the star on the log L - log T^ff curves as if the evolution were three times as 
slow as H-powered evolution. This provides a fairly good approximation of the behavior of the 
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models constructed by Iben (1984). The fraction of stars experiencing a final flash is not strictly 
determined by theory. Iben has calculated the probability that a star ignites helium in various 
post-AGB phases (Table 2, Iben 1984), obtaining different guesses for FF occurrence in 10 to 21 % 
of all stars in this mass range, and up to about 60 % when he includes stars that leave the AGB 
burning helium. Our parametrization is comparable with Iben's prediction. In fact, if we were to 
chose that 60 percent of all stars (in this mass range) have the chance to experience a final fiash, 
then about 12 % of the PNNi are in He-burning phase as observed in a synthetic HR diagram. 



2.6. The duration of the planetary nebula phase 

Following common wisdom, the duration of the FN phase is ^ 30,000 yr (e.g., Phillips 1989, 
and references therein). This is derived from the size of the largest observed PNs ('-^ 0.7 pc) 
coupled with a typical nebular expansion velocity of ~25 km s^^. In practice this assumes that all 
PNs remain visible as such until the expanding nebula has reached the maximum observed size, or, 
equivalently, that all PNs evolve in nearly the same way. On the other hand, the mere fact that PNs 
are produced by precursor stars in such a very extended range of initial masses (0.85 < Mi < 9 Mq) 
makes most unlikely that all PNs have the same lifetime. Simple arguments suggest that the 
time imax after the cessation of the superwind during which a FN remains observable scales as 
tmax oc ^^m\n^ ^ where MpN is the nebular mass and S'i^min is the minimum surface brightness 
for a FN to be detected (F2,F5). In the adopted mass loss parameterization, MpN corresponds to 
the mass ejected during the AGB superwind phase, and ranges from ~ 0.02 M0 for Mi = 0.85 M©, 
to over 1 Mq for Mi = 8 Mq . This factor of ~ 50 in FN mass therefore translates into at least a 
factor ~ 5 in tmax- To get the actual duration of the FN phase we must subtract to ^max the AGB 
to FN transition time, during which the object is a proto-FN, i.e., tp^ — ^max — ^tr • If '^tr ^ ^max? 
then when the central star has reached 10,000 K the nebular material has already dispersed, no 
observable FN is produced, and one has a so-called lazy post-AGB remnant (P2). To explore the 
effects of this subtle interplay between central star and (partially decoupled) nebular evolution, we 
have explored the case of mass dependent tmax, scaled with MpN- The nebular mass was derived by 
imposing its dependence on the plateau luminosity, since it is at the reach of a critical luminosity 
that the shell is ejected. The parametrization of MpN and tmax arc described in the Appendix (Eqs. 
Al and A2). Simulations of post-AGB populations with constant maximum FN ages will be also 
shown in this paper. 



2.7. The Montecarlo procedure 

Our population synthesis code starts with the option for a time limited or luminosity limited 
sample. The former option is to explore post-AGB populations within a fixed time range, (we are 
talking here about the stellar evolutionary time, not the duration of the PN phase, described in 
§2.6). This option is used to deduce the mass distribution of the PNN as derived from the synthetic 
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population analysis. The latter option, the luminosity-limited sample, is used to compare directly 
synthetic and observed diagrams. All stars in the synthetic population would have L> Lp (we 
assume Lf=1.0 Lq)^. 

To build the synthetic population, we proceed as follows: 

2.7.1. Time limited sample 

1) First, a value of the initial mass Mi is extracted, following the distribution: 

V'(Mi) oc Mi-(^+^) (13) 

2) The extracted value of M[ is entered into the initial mass-final mass relation to get the mass 
M of the post-AGB object. 

3) To proceed, we need to evaluate (see §2.4), and correspondingly the wind time scale 
can be calculated (see Eq. 9). 

4) A random time is then extracted, within a chosen range (0 < t < iiim)- 

5) If the extracted time t is larger than the wind timescale, we put t = t — t^, and this value 
of t is entered into the routines of §2.1 thus getting the corresponding location in the HR diagram. 
If t is smaller than the star is still in the wind phase, and we assign Te=3.7 and L = LpL. We 
call the resulting simulated star a wind object^. 

If we do not consider a FF, we go to point (8). 

6) Should we consider the effects of a final helium flash, another value of time tpF is extracted, 
with < tpF ^ himy so as to get a random value for the time at which the flash takes place 

7) If t < tpF) then this t is entered into the routines of §2.1 to get the luminosity and tem- 
perature of the star (in this case we have extracted a post-AGB star which will experience a final 
flash, but the flash has not taken place yet). 



Naturally, the luminosity limit can be varied accordingly to the observed situation we may want to reproduce. 

® Given our assumption on the MLR, we distinguish between wind objects and proto-PNN (in [9], below): wind 

objects arc those transition objects whose evolution is wind dominated, proto-PNNi arc all the other transition objects, 
located beyond the AGB. The assigned temperature for the wind object is rather arbitrary, but homogeneous to our 
mass loss choices. Even if we were to locate the wind objects within the transition area, the results would be not 

majorly affected, given the paucity of the samples. In a future paper, different mass loss choices will analyze the 
possible differences. Observationally, it may be hard to distinguish between these two types of objects. 

^°Note that Iben (1984) has calculated slightly different probabilities for FF occurrence at different phases, but we 
use a flat probability for simplicity. 
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Ift> tpF, then the time elapsed after the the flash {t' = t — tpp) is entered into Equation (12) 
to get luminosity of the post-FF star. To calculate the temperature we use the routines of §2.1. 

8) The transition time ttr is obtained following the prescriptions in §2.4. 

9) The object is finally classified into one of the following classes, and plotted with a different 
symbol for each class. The object is classified as a proto-Planetary Nebula Nucleus (proto-PNN) 

if t < ttr • It is classified PNN if t + ttr < tniaxj or t' + tpp + ttr ^ ^maxj with the further 
distinction between objects that have experienced a final flash (PNNHc) and those which did not 
(PNNH). It is classified as a post-PNN object otherwise, again distinguishing post FF and no post 
FF objects. 

2.7.2. Luminosity limited sample 

To build a luminosity-limited sample, we follow a similar procedure as in the time-limited 
sample, except in (4), the extraction of a random time, that is chosen in the range < t < 
tF{M) + 1^, where tp is the fading time to the minimum luminosity, L = 1.0 L© (see A3). 

The we proceed as in points (5) through (8) of the procedure in § 2.7.1, with the difference 
that tpF is extracted within the new time limits. 

3. Results 

Our code is flexible, suitable to produce post-AGB population synthesis for many applications 

and many sets of parameters. Far from being exhaustive, this paper includes only a small part of 
the possible applications. We have explored the issues and questions that we though to be of wide 
interests, and allow to advance in out post-AGB evolution knowledge. Other applications will be 
developed in the future. 

3.1. Synthetic diagrams and mass histograms: the basic model 

The first synthetic population shown here is a luminosity-limited sample of post-AGB stars. We 
extract 1500 objects, and we separate the wind objects, the PNNi, the proto-PNNi, and the post- 
PNNi following the prescriptions of Chapter 2. In Figure 15 we show the synthetic HR diagram, 
where 313 objects of the 1500 extracted are in the PNN phase, 21 are in the wind phase, and 9 
are proto-PNNi. The remaining objects are, following our prescription, already in the post-PNN 
phase. 

In the simulation of Figure 15 the residual envelope mass, M^, is chosen to be function of the 
core mass, as described in §2.4; we have used Weidemann's IMFMR, Salpeter's IMF, and no FF 
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objects (i.e., all PNN are hydrogen burners). The maximum PN time has been set following Eq. 
A2, with K=4xl0^, a relatively long time for PNe to disappear. We use this time through §3.5, 
considering furthcrly the effects of the maximum PN age in §3.6. We define the simulation plotted 
in Figure 15 as basic, and we will explore the parameter space by changing one parameter at the 
time with respect to this basic population. 

As expected, the synthetic stars cluster toward the lower masses, only the very low luminosity 
part being populated by higher masses. This is exactly the effect we see in all complete galactic 
PNNi samples. 

For more insight, let us plot the corresponding distribution of the visual magnitudes versus evo- 
lutionary timescale (Fig. 16). We translate temperatures and luminosities into visual magnitudes 
by following the bolometric correction from Code et al. (1976). As already evident in the previous 
plot, there are two mainly populated loci of the My — t plane as well, corresponding one to plateau 
luminosity of the low-mass objects, and the other to the stellar crowding toward low luminosities 
of the intermediate- and high-mass objects. The large fraction of PNNi that have decayed to post- 
PNNi at late evolutionary times is very evident in this Figure. This effect is mass-dependent, and 
smaller mass PNN disappear earlier in their life. 

In Figure 17 we show the (logarithmic) mass distribution for the basic synthetic population. 
This distribution has been obtained by running a sample of 1500 post-AGB stars with maximum 
HR diagram-life time of 30,000 years. The logarithmic scale has been chosen for a better view 
of the distribution in the whole mass range. The mass distribution of Figure 17 shows the clear 
clustering of model post-AGB stars around ~0.6 Mq, and a gradual spread to higher masses. 

Due to the particular choice of the time interval explored (t < 30,000 yr), we do not obtain 
post-PNNi with this time-limited selection. In Table 2 we show the composition of a selection 
of the synthetic populations shown in this paper. The columns list, respectively, the number of 
PNNi, wind objects, post-PNNi, and proto-PNNi stars for each extraction. Composition of typical 
populations may vary from one random extraction to another. From Table 2 we can directly 
compare the time-limited to the luminosity-limited samples. 

3.2. Effects of the initial mass function 

To test the variations of post-AGB distribution for different Initial Mass Function we have 
run synthetic populations with the following IMFs (see §2.2 for meaning of variable x and for SF 
choice): (1) the standard, Salpeter's (1955) IMF, with (1 -|- x) = 2.35 in the whole mass range 
considered; (2) Miller & Scab's (1979) IMF, with {1 + x) = 1.4 for M < IMq and (1 x) = 2.5 
for M > 1 Mq; (3) the IMF by Kroupa et al. (1991), with {I + x) = 1.85 for M < 1 Mq (this 
distribution is not defined for larger masses), and (4) an updated empirical parametrization by 
Scab (1998) with (1 + x) = 1.2 for M < IMq and (1 + x) = 2.7 for M > IMq. 
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In Figure 18 we show HR distributions accounting for IMF variations, when the IMF is the 
only parameters changing across the panels. In each case, we have extracted a luminosity-limited 
sample of 1500 post-AGB stars, with Weidemann's IMFMR, and mass-dependent residual envelope 
mass. The stars do not experience a final helium shell flash. In the case of Kroupa et al. 's low 
mass IMF, we use Scalo's (1998) IMF for masses larger than solar. 

By examining Figure 18 we see very little changes on the log L - log Tefj plane for different 
IMF in the mass range considered. The PNN populations in the case of variable IMF index (Fig. 
ISbcd) are larger for low masses respect to the basic model. The proto-PNN population changes 
slightly (we are dealing with low number statistics, anyway), but overall the distributions are 
similar in all four cases of Figure 18. Better to say, that the HR diagram is not the ideal way 
to pick up these differences, occurring in this case at low-luminosities, in a crowded part of the 
diagram. Observations of real situations such as those of the four panels of Figure 18 would not be 
distinguishable from one another. 

In Figure 19 we illustrate the mass distribution for the synthetic populations corresponding to 
Figure 18, each panel showing a different IMF. Naturally, time-limited samples have been used. It 
is worth noting that peaks in the distribution of less than a can be produced by the randomness 
of the simulation. 

The effects of the IMF are noticeable given the large sample of model stars. For example, 
Scalo's IMF (Fig. 19d) would produce a broader distribution toward the higher masses. Note that 
the logarithmic scale partially hides the fact that all these mass distributions are extremely peaked 
around 0.55 Mq. 

3.3. Effects of the initial mass-final mass relation 

As introduced in §2.3, we will use four different IMFMR and test their effects on the post-AGB 
populations. In Figures 20 and 21 we plot the results relative to (a) Weidemann's, (b) Ciotti et al. , 
(c) Herwig's, and (d) IR's IMFMRs. Figure 20 shows the distributions on the HR diagram. On 
this plane, the sample with IR's IMFMR stands aside, since the IR's prescription produces more 
massive post-AGB. 

The differences among the other distributions are more evident in Figure 21, where we plot 
the visual magnitude versus evolutionary time of the synthetic stars. By comparing the basic and 
Herwig's IMFMR diagrams, above Mv=5 the two distributions are similar, but the differences 
appear and get more evident for fainter objects. 

To compare these results with the actual observed stars is beyond the purpose of this paper. 



^^In Fig. 18, and after, we use the following conventional nomenclature for figure panels: a=lower left, b=upper 
left, c=lower right, d=upper right. Panel (a) always shows the basic sample. 
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Nonetheless, we can state that the different IMFMRs could be inferred mostly at low brightness, 
with consequent higher uncertainty of the comparison with the data. 

The mass distributions of Figures 22 set aside once again the IR IMFMR. This old parametriza- 
tion allows a continuum of masses between 0.55 and 1.4 Mq. The other 3 cases have similar low- mass 
distributions, while the high-mass distributions are similar in the Ciotti et al. 's and the Herwig's 
IMFMRs, as the basic sample allows for few very high mass stars. The observational consequences 
of the four different IMFMRs are mainly constrained in the high-mass tail of the distributions. 

3.4. The residual envelope mass 

The envelope mass that is left on the star after the envelope ejection, Mj^, plays a fundamental 
role in the following stellar fate. Our basic model uses M^=f(M), with an ad hoc parametrization 
(see §2.4). Since the value of Mf- is undetermined, as well as its distribution with respect to the 
other physical parameters, we chose to explore scenarios of post-AGB evolution both with with 
variable and constant M^. We also produce populations with random residual envelope masses, 
to represent the extreme (but not unrealistic) indetermination of M^, and consequently of the 
transition time. 

In Figure 23 we show the results with the different scenarios: the basic model (a) has mass- 
dependent M^, model (b) has random (with M^<0.1 Mq), model (c) also has random 
Mf- (with M^< 1 X 10""^ Mq), and model (d) has constant Mf-=5 x 10"^ Mq, independent on 
the progenitor mass (lazy-AGB stars have not been produced in these simulations, given the choice 
of K in equation A2) . 

The effects of the indetermination of the residual envelope mass are striking on the resulting 
post-AGB populations. A complete description on mass loss for low- and intermediate-mass stars 
does not exist to date. Our experiments (e.g., Fig. 23) represent a way to determine the effects of 
the mass dependence on the MLR on observable sets of PNNi and related objects. Once again, in 
this paper we do not compare these effects with the observed stars, we just set the stage for future 
comparisons. When using M^=5xlO~^M0, the distribution is very similar to the basic case. 

We use the random M^ to further analyze the effects of the mass loss on observed populations. 
The distributions that we obtain as a result of the simulations with random Mf- are very different 
depending on the upper limits that we set for M^. For instance, for M^< 1 x 10~^ Mq, no wind 
nor proto-PN are created. Furthermore, very few high-luminosity PNe are available with this low- 
limit random M^ (see also Table 2). These effects could be observed in PNNi and proto-PNNi at 
known distances. At lower luminosities, the two distributions look similar especially if we take into 
account observing errors. 

The mass distributions of post-AGB populations with different choices of M^ are shown in 
Figure 24. The effects of changing M^ are even more extreme in these time-limited samples. 
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Let us examine Figures 24b and 24c (upper- left and lower-right panels). In both panels, Mf- has 
random values, the difference being the maximum allowed M^. In the simulation of 24b, where 
the maximum residual envelope mass is 0.1 Mq, the distributions is dominated by the objects in 
the wind phase, and the lack of PNNi (actually, the simulation produces two PNNi, hiddent in the 
log-representation). The most notable characteristics of panel 24c are the lack of wind objects or 
proto-PNNi. The very low Mf- make the transition time very short, and, as a consequence, the 
lack of proto-PNNi. Studies of a sample including OH/IR and other transition objects, and PNNi 
within the same environment, should reflect the mean value of in the sample. 

In Figure 24d most post-AGB stars belong to the pre-Planetary Nebula (wind and proto-PNN) 
phases. Only a few percent of the stars are PNNi. The situation of Figure 24d has been produced 
only to show an extreme case, since it is highly unrealistic that all stars have the same residual 
envelope mass, independent of their initial mass on the main sequence, and on their core mass. 

3.5. The final helium-shell flash 

In this paper we did not use the helium-burning (nor a combination of hydrogen and he- 
lium burning) stellar models to determine the effect of a helium-burning PNN population. The 
parametrization that we use is described in §2.5, and it agrees with observations. Following, we 
show how the presence of helium-burning stars can effect the observable post-AGB populations. In 
Figure 25 we show the synthetic population on the HR diagram for a sample of 1500 post-AGB 
stars in which 20 % of the post-AGB stars have experienced a final flash. The fraction of stars in 
post-FF is chosen in agreement with the guesses based on the evolutionary models by Iben (1984). 
We compare this population with the basic model. We note, from Table 2, that the overall compo- 
sition of the sample in terms of ratio of PNNi to other components is (statistically) very similar to 
the basic sample. 

In Figure 25 we have separated the He- and H-burning stars. Panel (a) shows the basic model 
(1500 objects with Salpctcr's IMF, Weidemann's IMFMR, Mf=i{M), and no post-FF). In the 
other panels, we show (b) the post-AGB stars having experienced a FF (20% of the sample); (c) 
stars that are in H-burning, post H-burning and/or will experience FF (the latter group of stars is 
observationally indistinguishable from H-burning stars); (d) the composite population of panels (b) 
and (c) , thus the observable post-AGB population. Figure 26 shows the same populations of Figure 
25 in the My — t plane. The simulated He-burning populations stand out for lack of low-brightness 
PNNi. 

It is not easy to observationally discern among H- and He-burning stars. The problem is that 
the stellar abundances are not easy to measure in PNNi. In general, H-depleted post-AGB stars 
are believed to be He-burners, but it is not so easy to single out the H-burning stars, as H-rich 
post-AGB could indeed be He-burners as well as H-burners. The simulation of Figures 25 and 26 
can help us in this respect. In fact, the ratio of PNNi brighter than, say, 5 magnitudes in the basic 
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sample is about a third of the total sample, while in the composite population (Fig. 26c) this ratio 
goes up to half the sample. A homogeneous, complete observed PNN sample should show this kind 
of discrepancy. 

3.6. The planetary nebula life time 

Through this paper we have assumed that the maximum age for a PN depends on the FN mass, 
as in Equation A2, with constant K=4xl0^. We have also produced other synthetic populations by 
keeping the nebular mass dependence, and changing the constant. If we keep all other parameters 
as in the basic model, we obtain that the number of PNNi goes to zero as K is lowered to about 
5x10^. This happens because the shorter lifetime of PNe makes most of the post-AGB population 
be in a post-PNNi phase. In this case, a large number of lazy AGB stars are also produced. In 
general, lazy AGB stars are produced by running the basic model with K< 8 x 10^. 

If we were to eliminate the nebular mass dependence from tmax the situation changes noticeably. 
In Figure 27 we show a simulated population with all parameters identical to the basic model, except 
we have fixed the PN life time to 30,000 years for all nebulae. The result is remarkably different 
than the basic model (see Fig. 15 and Table 2). Most low-mass post-AGB stars are in post-PNN 
phase. If we were to lower the maximum PN life time to 10,000 years (not an unreasonable choice) 
we obtain that most PNNi in the diagram would have a high mass. This is in contrast with the 
observations, thus our simulations lead us to believe that the lifetime of PNe depends on the nebular 
mass. Detailed comparison with homogeneous dataset should be used to confirm this inclination. 

4. Summary and future work 

We have used up-to-date evolutionary tracks as templates to build a code for post-AGB popula- 
tion synthesis. Our models aim at understanding the fine tuning of post-AGB evolution, including 
the consequences of IMF and IMFMR, the transition time and its correlation with the residual 
envelope mass, the actual duration of the PN phase, the relative population of proto-PNNi, wind 
objects, and post-PNNi stars, and the occurrence and timing of the FF phase. Our synthetic 
tracks, available to obtain log Tpg and log L for post-AGB stars of any stellar mass in the range 
0.85 < Mi < 9Mq, arc obtained with very high precision in reproducing the actual evolutionary 
tracks. The interpolation method is simple and fully explained in this paper, so that the readers 
can reproduce the synthetic populations included here. 

In this paper we have shown a sample of the possible applications, without going into detailed 
comparisons with the observed data. Among the results shown here, we found that (1) the de- 
pendence of the synthetic populations on the assumed IMF and IMFMR is mild; the post-AGB 
populations are not ideal indicators of the IMF. (2) The residual envelope mass, after the envelope 
ejection, has a strong effect in determining the subsequent post-AGB evolution; its indetermination 
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produces very high indetermination on the transition time, and ultimately on the resulting post- 
AGB populations. (3) The ratio of He- to H-burning PNNi can be reproduced with population 
synthesis. 

The central importance of this paper consists in showing that the many fundamental variables 
of stellar evolution have a major role in determining post-AGB populations. The variation and 
indetermination of these parameters should not be overlooked in comparing data and theory. 

The theoretical work contained in this paper was implemented with the goal of being versatile 
and useful for a full host of applications. Among other possible uses of these synthetic tracks 
are the studies of evolutionary effects on the Planetary Nebula Luminosity Function (PNLF), and 
how those translate in the variation of the extragalactic distance scale, as derived from the PNLF. 
Stanghellini (1995) has shown in a preparatory work that the PNLF is at variation with respect 
to the transition time. The updated models presented here allow the most detailed study of the 
PNLF on the evolutionary parameters. 

Future applications also include the simulations of bulge, elliptical galaxy, and Magellanic 
Cloud populations, where a different treatment of the SF/IMF should be used, and possibly different 
stellar chemistry. 

Thanks to Laura Grcggio for her help in the interpolation techniques, to Antonclla Nota and 
many others for their bibliographic indications, and to an anonymous referee for pointing out an 
error in the previous draft and for many suggestions. 



-20- 



A. Appendix 

In order to parametrize the duration of the PN phase, we assume that the nebular mass is 0.02 
Mq for M=0.55 Mq (core mass), and 1.5 Mq for M=1.2 Mq. The planetary nebula mass can 
be described by the following correlation between nebular mass and plateau luminosity: 

MpN = 4.495 X 10-^LpL^-^^^, (^1) 

and 

W = i^MpN (^2) 

To calculate the fading time to L = I.OLq, tp, we use the VW tracks and we interpolate versus 
the (post-AGB) mass, obtaining 



tF{M) = 7.0957 - 2.1441M. 



(^3) 
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Fig. 1. — Effective temperatures versus evolutionary time for M=0.569, 0.597, 0.633, 0.677, 0.754, 
and 0.9 Mq. 

Fig. 2.— Lumonisities versus evolutionary time for M=0.569, 0.597, 0.633, 0.677, 0.754, and 0.9 

Fig. 3. — Effective temperature versus normalized age in the hydrogen burning phase. 

Fig. 4. — Effective temperature versus normalized age in the quenching phase. 

Fig. 5. — Normalized effective temperature versus normalized age in the cooling phase. 

Fig. 6. — Normalized effective temperature versus normalized age in the white dwarf phase. 

Fig. 7. — Luminosity versus normalized age in the Hydrogen burning phase (or plateau) . 

Fig. 8. — Normalized luminosity versus normalized age in the cooling phase, for M=0.569, 0.597, 
0.644, and 0.644 Mq. 

Fig. 9. — Luminosity versus normalized age in the white dwarf phase. 

Fig. 10.— The IMFMR: from IR (solid hne), Ciotti et al. (1991, long dashed line), Weidemann 
(1987, short-dashed line), and Herwig (1997, dash-dotted line). 

Fig. 11. — The transition time (thick solid line), obtained with as a function of the core mass, 
as explained in §3.1. Also plotted are the wind (thin solid line), and nuclear (dashed line) time 
scales. 

Fig. 12.— As in Figure 11, but with Mg^=3xl0-^ Mq. 
Fig. 13.— As in Figure 11, but with Me^=5xl0~^ Mq. 
Fig. 14. — As in Figure 11, but with random < 5e — 3 Mq. 

Fig. 15. — The basic post-AGB synthesis, luminosity- limited. Solid circles = PNNi; squares = wind 
object; open circles = proto-PNNi. The tiny dots show the locations of the post-PNNi. Tracks 
correspond to M=0.535, 0.569, 0.597, 0.633, 0.677, 0.754, and 0.9 Mq. 

Fig. 16. — The basic post-AGB synthesis (luminosity-limited distribution), on the My — t dia- 
gram. Solid circles = PNNi; squares = wind object; open circles = proto-PNNi. The tiny dots 
show the locations of the post-PNNi. Tracks correspond to M=0.535, 0.569 and 0.9 Mq. 

Fig. 17. — Logarithmic mass distribution for the basic sample (time-limited distribution). Solid 
line = PNNi; broken line = proto-PNNi; shaded histogram = wind objects. 

Fig. 18. — Synthetic post-AGB populations on the HR diagram. The IMFs are from (a) Salpeter 
(1955) (b) Miller & Scalo (1979); (c) Kroupa ct al. (1991); and (d) Scalo (1998). Other parameters 
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are as in the basic model, which is panel (a). Symbols and tracks as in Fig. 15. 

Fig. 19. — Mass distribution of synthetic populations with different IMFs. The IMFs are from (a) 
Salpeter (1955); (b) Miller & Scalo (1979); (c) Kroupa et al. (1991); and (d) Scab (1998). Other 
parameters are as in the basic model. Histogram lines as in Fig. 17. 

Fig. 20. — The log T^g -log L plane population synthesis for (a) Weidemann's, (b) Ciotti et al. 's, 
(c) Herwig's, and (d) IR's IMFMR. All other parameters are as in the basic models. Symbols and 
tracks as in Fig. 15. 

Fig. 21.— The My - t population synthesis for (a) Weidemann's (1987), (b) Ciotti et al. 's (1991), 
(c) Herwig's (1997), and (d) Renzini & Voh's (1981) IMFMR. All other parameters are as in the 
basic model. Symbols and tracks as in Fig. 16. 

Fig. 22. — The mass distributions of populations with different IMFMR: (a) Weideman's, (b) Ciotti 
et al. 's, (c) Herwig's, and (d) Renzini &; Voli's IMFMR. All other parameters are as in the basic 
model. Histogram lines as in Fig. 17. 

Fig. 23.— The My-t relation for different choices of the M^; (a): M^=i{M); (b): random M^, 
with Mg^<0.1 Mq; (c): random M^, with M^< 1 x 10~^ M©; (d): constant Mg^=5xl0~^ Mq. 
All other parameters are as in the basic model. Symbols and tracks as in Fig. 16. 

Fig. 24. — Mass distribution for the choices of Mj^ as in Fig. 23: (a): M^=f{My, (b): random 
M^, with Me^<0.1 Mq; (c): random M^, with M^< 1 x 10"^ Mq; (d): constant Me^=5xl0-^ 
Mq. All other parameters are as in the basic model. Histogram lines as in Fig. 17. 

Fig. 25. — The effects of a final helium-shell flash during cooling, on the HR diagram, (a) basic 
population; (b) post-FF stars; (c) H-burning or post H-burning stars; (d) population of post-AGB 
stars with 20% objects in post-FF phase (i.e., b-l-c populations). All other parameters are as in the 
basic model. Symbols and tracks as in Fig. 15. 

Fig. 26. — The effects of a final helium-shell flash during cooling, in the My — t diagram, (a) basic 
population; (b) post-FF stars; (c) H-burning or post H-burning stars; (d) population of post-AGB 
stars with 20% objects in post-FF phase (i.e., b-|-c populations). All other parameters are as in the 
basic model. Symbols and tracks as in Fig. 16. 

Fig. 27. — The post-AGB synthesis, luminosity-limited, with PN mass independent maximum PN 
lifetime (tmax=30,000 yr). Symbols and tracks as in Fig. 15. 
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Table 1: Characteristics of evolutionary tracks phases 



Phase 


M 




log Tgff eni-i 


log Lend 




[Mq] 


[yr] 


'[K] 


L W J 


H-huminq 


0.569 


3.215x10^ 


4.907 


1.47 




0.597 


1.098x10^ 


5.153 


3.19 




0.633 


4448 


5.238 


3.348 




0.677 


2667 


5.292 


3.474 




0.754 


1067 


5.39 


3.699 




0.9 


141.7 


5.581 


4.02 


quenching 


0.569 


3.614x10^ 


5.038 


2.579 




0.597 


1.214x10^ 


5.071 


2.458 




0.633 


4862 


5.112 


2.424 




0.677 


2925 


5.162 


2.527 




0.754 


1177 


5.315 


3.033 




0.9 


167.1 


5.403 


2.953 


cooling 


0.569 


3.031x10^ 


4.907 


1.47 




0.597 


3.258x10^ 


4.9 


1.297 




0.633 


3.113x10^ 


4.906 


1.203 




0.677 


3.165x10^ 


4.911 


1.125 




0.754 


3.318x10^ 


4.91 


0.919 




0.9 


2.21x10^ 


4.952 


0.855 


white dwarf 


0.569 


3.068x10^ 


4.659 


0.112 




0.597 


4.044x10^ 


4.62 


-0.127 




0.633 


4.681x10^ 


4.596 


-0.299 




0.677 


4.973x10^ 


4.583 


-0.411 




0.754 


7.012x10^ 


4.547 


-0.683 




0.9 


1.164x10^ 


4.545 


-0.877 
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Table 2: Composition of typical synthetic post-AGB populations 



run 


rJNJNi 


Wind objs. 


post-FJNJNi 


proto-rJNJNi 


Luminosity-limited sample {L > ILq) 










basic 


313 


21 


1157 


9 


TV yr*n o 1 

Miller & Scalo 


246 


18 


1230 


6 


Kroupa et al. (1991) 


265 


13 


1212 


10 


bcaio (1998) 


z8z 


14 


1199 


5 


Ciotti er m. (1991) 


401 


10 


1083 


D 


TT " / 1 i \( W \ 

Herwig (1995) 


352 


12 


1125 


11 


Renzini & Voli 




1 K 


ovo 





random Mf{< 0.1 M©) 


191 


520 


784 


5 


random Mf-{< 1 x 10"^ Mq) 


226 




1274 





M^^= 5 X 10-3 Mq 


232 


41 


1225 


2 


20 % post-FF 


275 


11 


1210 


4 


*max=30,000 yr 


52 


17 


1426 


5 


Time-limited sample {tum < 30, 000 yr) 










basic 


940 


397 




163 


Miller & Scalo 


959 


400 




141 


Ciotti et al. (1991) 


975 


381 




144 


Herwig (1995) 


970 


288 




142 


random M^{< 0.1 M©) 


2 


1498 






random M^{< 1 x 10"^ Mq) 


1215 




285 




M^=5 X 10-=^ Mq 


578 


813 


8 


101 
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